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On the basis of a microscopic theory, the signatures of many-particle correlations in Two- 
Dimensional Fourier- Transform Spectra (2D-FTS) of semiconductor nanostructures are identified 
and compared to experimental data. Spectra in the photon energy range of the heavy-hole and 
light-hole excitonic resonances show characteristic features due to correlations, which depend on the 
relative polarization directions of the excitation pulses. 

PACS numbers: 78,47.-|-p, 42.50.Md, 42.65. Re, 78.67.-n 

Recent reports illustrate the potential of a novel method, known as "Two-Dimensional Fourier- Transform Spec- 
troscopy" (2D-FTS), for the investigation of many-particle induced correlations in semiconductor structures 0, [2, S]- 
2D-FTS is based on a four- wave- mixing experiment, where three excitation pulses are separated in time and is hetero- 
dyne detected to fully characterize its phase. The signal is transformed into frequency domain wt (emission frequency) 
and LOr (excitation frequency), both with respect to real time t and time separation r of the first two pulses, respec- 
tively. The third pulse is delayed with respect to the second one by T. 2D-FTS is widely used to study vibrational 
0, i, i and electronic excitations 0, i i in molecules. Applying 2D-FTS in the optical regime we investigate the 
complex interplay between exciton, biexciton and continuum excitations, in semiconductor nanostructures. 

In previous publications demonstrating 2D-FTS of semiconductors 0, 0, 13| , the experimental results were compared 
to a phenomenological theory based on extending the Optical Bloch Equations to include terms that describe excitation 
induced dephasing (EID) and excitation induced shift (EIS) [12]. Here, we use a microscopic theory to calculate 

2D-FTS. Furthermore, the model predicts that the 2D-FTS qualitatively depend on the polarizations directions of 
the excitation pulses. Our results demonstrate that for 2D-FTS, the influence of many-particle correlations on the 
spectra can clearly be identified and are in agreement with the experimental findings. 

For the qualitative modeling of 2D-FTS in such systems and in order to keep the numerical requirement within 
reasonable limits we use a microscopic many-body theory and apply it to a one-dimensional tight-binding model. 
Therefore, a quantitative agreement between theory and experiment can not be expected. It has, however, been shown 
that many important signatures of nonlinear optical experiments performed on quantum wells can qualitatively well 
be reproduced by this model [H, [H, [H, [H O [H . 

The optical nonlinearities are treated up to third order in the coherent x^^-'-limit beyond the Hartree-Fock level. 
In order to separate the correlation effects from those due to the first-order Coulomb interaction (i.e., Hartree-Fock) 
we use a set of equations for the interband coherence P and two-exciton amplitude B which in symbolic form reads 

dP 

- ih— = h{ujx - il)P + Vc P* PP - VcP*B + fi*E 



dt 



dB 



- ^i*E P* P, (1) 



-ih— = h{LUBx+iP)B + PPVc. (2) 

In fact, all terms in the microscopic version of the equations depend on the spatial indices of the tight-binding model 
as well as band indices, referring to electron, heavy-hole (h) and light-hole (1) bands, respectively. Details of the 
model and the set of microscopic equations that are used in our numerical evaluations are presented and discussed in 
Ref. [l^. In Eqs. ^ and ([2]) Vc is the Coulomb matrix that is responsible for the existence of exciton and biexciton 
resonances and interactions between them. In Eq. ([1]) for the polarization P, the Hartree-Fock terms are given by the 
kinetic term hux , the renormalization term Vc P* PP and the coupling to the light field containing the Pauli-blocking 
IJ*E P* P and the linear term iJ*E. The correlation term beyond Hartree-Fock is given by VcP*B. As given in 
Eq. ([2]) the biexciton correlation B is driven by PPVc- The biexciton energy is Hlobx- Coherent x^^"* -equations with 
h- and 1-transitions have been previously analysed, see e.g., [isl. [2^. 

The quantum well samples used in this study show heavy hole and light hole excitonic resonances. Their spectral 
positions are determined by the kinetic parts of Eqs. ([T]) and which also contain the phenomenological dephasing 
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rates, 7 = \/T2^'^'' for excitons and (3 = for biexcitons. By such rates we can distinguish biexcitons 

formed by h-excitons {hh), by 1-excitons {II) and mixed h- and 1-excitons {hl,lh). The optical transitions are given 
by dipole matrix elements ^ chosen to model the relevant selection rules for dipole transitions in III-V semiconductor 
structures [H, This model agrees with the selection rules for the nonlinear polarization presented by M. Lindberg 
et al. for the three beam situation [2^ . 

Based on this model we present microscopic calculations for 2D-FTS for semiconductor structures. We choose 
realistic material parameters, ratio between h- and 1-effective masses and dipole matrix elements, energetic h- and 
1-offset, and Coulomb strength. The experimental conditions such as pulse duration, detuning, dephasing times, 
polarization directions and the sequence of the pulses are incorporated into the model. The dephasing times of the 
various excitons and biexcitons are taken as phenomenological parameters. 

In this paper, we concentrate on the situation where the first pulse enters the nonlinear polarization conjugated, 
the so called rephasing case Q. The figures show the amplitude \E{uJt,i^T,T)\, where the electric field E{u!t,i^T,T) 
is proportional to the third-order polarization multiplied by the imaginary unit, i.e., iP{uJt,WT,T). The emission 
frequency, ujt, is given on the horizontal-axis. The emission frequency is used to define the sign of frequency, thus 
the conjugation of the first pulse means that the excitation frequency, cUr, is negative. It is given on the vertical axis. 
The h- and 1-exciton peaks then show up on the diagonal in the left-upper (low energy) and right-lower (high energy) 
corners, respectively, while the off-diagonal peaks are due to coupling between h- and 1-exciton resonances. 

For the particular system considered in this paper, i.e., a GaAs quantum well, the nonlinear optical properties 
have been explored quite extensively in the past using four-wave-mixing and pump-probe experiments. This system 
therefore presents a good testing ground for the exploration of the features of 2D-FTS. In particular, the sample 
used in the experiment is a 10 period multiple quantum well structure with 10 nm GaAs wells and 10 nm AlGaAs 
barriers. It is held at 8K. Experimental spectra are obtained using the apparatus described by Zhang et al. with 
the addition of polarizing optics. 

In the following we concentrate on the amplitude features of 2D-FTS. The distribution of the heights of the peaks 
strongly depends on the interplay of the material and experimental parameters, including tuning of the excitation 
pulses, shape and temporal full width at half maximum (FWHM) of the pulses, ratio iJ.^/fJ' of h- vs. 1-dipole matrix 
elements (which due to band-mixing in the quantum well structure depart from their bulk values) and dephasing 
times. We found that ~ 2.1 is a good choice for modeling the experiment. The simplest approach to model 

the h-l-exciton system is a three level system (V-system) without any interactions, but with coupling of two excited 
levels (modeling the h- and 1-excitonic resonances) to a common ground state. In this model it is not possible to 
obtain a peak height distribution in the 2D-FTS features that is asymmetrical with respect to the diagonal, if the 
excitation pulses are extremely short and fi^ = fiK The asymmetry of such spectra in our experiments on quantum 
well systems appears because of the interplay of several effects: different dipole matrix elements, dephasing times 
of excitons and biexcitons, degree of overlap with the spectrum of the excitation pulses and many-particle couplings 
within the system, as supported by the theoretical spectra. In our model the electric field of the excitation pulses is 
proportional to exp{—P/5^), where S is the Gaussian width of around 115 fs, corresponding to intensity FWHM of 
135 fs and a spectral FWHM of 13.5 meV. 

In Figure [1] we compare experimental (first column) and theoretical results for the amplitudes from the full calcula- 
tion (second column) and the Hartree-Fock part (third column) for cr+ (t+(t+, co- linear XXX, and cross-linear YXX 
polarizations. In all figures the central excitation energy is ImeV above the 1-exciton in order to compensate for the 
small dipole matrix element of the 1-exciton. 

Numerical results for the full calculation (second column) show agreement with experiments, not only for the 
h-, 1-exciton and mixed peaks, but also for the h-continuum at higher ujr (the vertical continuum contributions in 
the experimental data are more apparent in a previous publication 0). Note that in the theoretical figures, the 
h-continuum, appearing close to the 1-excitonic peaks as dominantly vertical structures, is decomposed into discrete 
peaks because only TV = 40 sites have been taken in the tight-binding approach for numerical reasons. It should be 
noted that the continuum appears dominantly as a vertical feature close to the mixed peaks in the lower-left corner. 
This demonstrates that the lower energetic part of the continuum can not be seen as a pure inhomogencous spectrum 
in the model containing the many-particle interaction (isl . [2^ . In an interaction- free model the continuum would 
in fact be purely inhomogencous, i.e., it would consist of independent two- level systems belonging to a spectrum of 
fc-vectors. If this were true, one would expect a feature of continua along the main diagonal close to the 1-excitonic 
peak in the right-lower corner. In the frame of the Hartree-Fock approximation, where the correlations are ignored, 
the features due to the h-continuum clearly show up at the higher energies for all polarization cases (Fig. [1] third 
column) . 

The co-linear situation (Fig. [T] middle row) looks similar to the co-circular case (Fig. [T] upper row) , as is to be 
expected from the selection rules. It is remarkable, that the upper right non-diagonal mixed peak for co-circular case 
(Fig. [TIj) is very small. As a comparison with a Hartree-Fock result (Fig. [TJ:) shows, this is not due to correlations, 
but results from an interplay of detuning, dipole matrix elements, and dephasing times. The mixed signature in 
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FIG. 1: 2D-FTS spectra amplitudes (normalized). Left column: experimental. Middle column: full x calculation. Right 
column: Hartree-Fock calculation. Upper row: a^a^a'^ polarization, T = lOOfs. Middle row: XXX co-linear polarization, 
T — lOOfs. Lower row: YXX cross-linear polarization, T — 350fs. Dephasing times are: — L3ps, — 0.8ps for excitons 
and T^'' = 1.3ps, Tj' = O.Sps, T^^ = T^'' = O.Sps for biexcitons. Theoretical spectra are weakly inhomogeneously broadened 
by Gaussian width of 0.3 meV. 

the left-lower corner for the Hartree-Fock calculation (Fig. (H:) shows discrete weak features due to the h-l-excitons 
and the h-continuum. The many-body correlations lead to the increase of this mixed signature at higher excitation 
energy and vertical features of the continuum as well. Unlike to the co-circular case, the influence of correlations 
on the mixed signature at higher excitation energy is less pronounced for the co-linear case. The Hartree-Fock 
approximation (Fig. [Tf) and the full calculation (Fig. [T^) do not differ that much for the co-linear case. This can be 
understood by consideration of a simple V-system, where for the excitation a~^a^cr^ the h- and 1-transitions resemble 
two independent two-level systems. The linear case contains both cr"'" and a~ transitions, thereby such excitation 
couples the h- and 1-exciton transitions resulting in two V-systems. This coupling generates the mixed signatures and 
appears for both linear polarized cases in the Hartree-Fock calculations (Fig. [If: XXX, i: YXX). 

The cross-linear polarized excitation case (YXX) shows additional signatures from two-exciton resonances (bound 
and unbound biexcitons), as predicted by the selection rules and supported by the theoretical results (Fig. [Th). This 
is due to a suppression of the excitonic features with respect to the biexcitonic ones JJJ. Boimd biexcitons show up 
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on the low emission-energy side and unbound ones on the high emission-energy side of the exciton [isl . 0, [T^ . The 
biexciton contribution results in a horizontal elongation of the peaks compared to the XXX-ca.se. The Hartree-Fock 
calculation (Fig. [TJ) clearly shows, that the horizontal elongation of the YXX spectra (Fig. [TJi) is due to correlations, 
i.e., bound and unbound two-exciton states. In the Hartree-Fock limit the horizontal elongation is absent (Fig. [TJ). 
Interestingly, the continuum contribution and respectively the elongation of the signatures at higher excitation energy 
is less developed as compared to the co-linear case, which is supported by the experiment (Fig. [TJi,g). For the cross- 
linear situation one can show ^231 that indeed the continuum becomes suppressed due to cancellation effects that 
result from Coulomb-correlations in this situation. Indeed, the Hartree-Fock calculation shows a more pronounced 
continuum contribution (and concomitantly of the superimposed 1-exciton peak) in the right-lower corner of the 
spectrum, as compared to the full x'^'^^-calculation. 

It is surprising that the lower left non-diagonal peak is the strongest one in the theoretical spectra for the YXX-case 
(Fig. [TJi), in contrast to the experimental data (Fig. [1^), where the h-peak is the strongest one. At present, the 
reason for this discrepancy is not clear. The peak in question is a mixture of three contributions, the h- and 1-excitonic 
resonances and the h-continuum. Therefore it is reasonable to assume an effective dephasing of this peak, which is 
not included in our present treatment. 

A comparison between the results of the full and the Hartree-Fock calculation shows that the left-lower non-diagonal 
peak is less prominent in the latter spectra. This difference suggests that correlations enhance this signature in the full 
treatment. These couplings could well be more sensitive to relaxation processes than the trivial ones which already 
appear in the V-system or within a model, which only considers the phase-space (Pauli-blocking) nonlinearities. 

We have not studied the dependence of the peak heights on the spectral shape of the incident laser pulses which we 
have here always taken as Gaussian. Our calculations have shown that the temporal duration and thus the spectral 
width of the laser beams strongly influences the 2D-FTS. Therefore, deviations of the pulse envelopes from a Gaussian 
shape will clearly significantly alter the peak heights. 

The amplitude of the non-diagonal signatures relative to the diagonal ones generally depends on T. Already for a 
V-system one finds that the two non-diagonal peaks show beating with respect to T, where the period is 2TTh/AEh,i, 
with AEh.i the energetic h- and 1-exciton separation. In the experimental data shown here T = 350fs has been 
taken. With AEhj — 11.5meV, corresponding to a beating period of 360 fs, this is close to the situation of maximum 
non-diagonal peaks. However, as mentioned above, the left-lower signature is definitely too strong if compared to the 
experiments. In our calculations for the YXX-ca.se we have used T = 350 fs. 

It should be noted that biexcitonic features can be observed in different kinds of nonlinear optical experiments for 
certain polarization cases. The 2D-FTS always show the presence of biexcitonic features in the plots for amplitude 
or real part (not shown here), albeit to different extent, depending on polarization. 

Since 2D-FTS is based on a photon-echo experiment, it is not surprising that it can be used to measure the dephasing 
time T2 even in the presence of inhomogeneous broadening. Specifically, inhomogeneous broadening will elongate peaks 
along the diagonal, while the width perpendicular to the diagonal is inversely proportional to T2. However, 2D-FTS 
is clearly superior to a traditional one-dimensional photon echo measurement when multiple resonances are involved, 
as is the case here. The photon-echo signal mixes the contributions of the resonances, possibly leading to ambiguity 
in assigning the dephasing times. This ambiguity does not occur in 2D-FTS as the peaks corresponding to distinct 
resonances are well separated. 

A traditional one-dimensional transient four-wave-mixing (or photon echo) measurement also has trouble yielding 
T2 when there are strong many-body correlations, which may result in a rapid decay of the signal unrelated to T2 
[2^. 2D-FTS is much better at revealing T2 of individual resonances. As an example we consider theoretical 2D FTS- 
amplitude plots for cross-linearly polarized pulses, which show the correlation-induced features most prominently. 
However, now the 1-transitions have been given the same dephasing time as the h-ones, i.e., = T2 — 1.3ps, 
while the theoretical results shown in Figure [T] correspond to T2 = 1.3ps and Tj = 0.8ps. Figure [5] shows, that the 
distribution of the peak heights changes due to longer dephasing time of 1-exciton. The comparison with Fig. [T}i 
demonstrates the large sensitivity of the 2D-spectra on T2-times. 

In summary, we have investigated correlation effects in quantum well systems by applying the Two-Dimensional 
Fourier- Transform Spectroscopy (2D FTS). This method provides a wide spectrum of information about the many- 
body correlations simultaneously, such as: the strength of the couplings between excitons, biexcitons and continua, 
character of continuum excitations, dephasing times, and dependences of these features on the polarization directions. 

A first comparison between experiment and theory has been performed on the basis of a microscopic model. 
Depending on polarization directions of the excitation pulses, characteristic signatures of many-particle correlations 
can be identified in the amplitude spectra. Differences between experiment and theory point towards the action of 
relaxation processes on many-particle correlations, which are not included in the present purely coherent treatment. 
Additional important information is contained in the real- and imaginary part of the spectra. Work is in progress to 
analyse these parts of the signal and to exploit these spectra for semiconductor structures in the energetic range of 
both excitonic and continuum transitions. 
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FIG. 2: Amplitude of 2D-FTS for YXX case. Dephasing times are: Ti, = = 1.3ps for excitons and Ti,' = T^^ = 1.3ps, 
rphi _ rpUl _ Q g5pg fQj. biexcitons. Theoretical spectra are slightly inhomogeneously broadened by Gaussian width of 0.3 meV. 
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